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We propose that a hierarchical spectrum of sterile neutrinos (eV, keV, 10 13-15 GeV) is considered 
to as the explanations for MiniBooNE and LSND oscillation anomalies, dark matter, and baryon 
asymmetry of the universe (BAU) respectively. The scenario can also realize the smallness of active 
neutrino masses by seesaw mechanism. 
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The compelling evidences from solar, atmospheric, re- 
actor, and accelerator neutrino experiments have estab- 
lished the phenomenon of neutrino oscillations. The 
standard description is that the experimental data can 
be nicely explained by the mixings between the flavor 
and mass eigenstates of the three neutrinos in Standard 
Model (SM), the so-called "active" neutrinos. The uni- 
tary mixing matrix is parametrized in terms of three ro- 
tation angles (#12, #23, #13) and one Dirac CP violat- 
ing phase 8c p- The probabilities of flavor oscillations 
are governed by the 6ij and two mass-squared differences 
Am 2 12 ~ 7.59xl(T 5 eV 2 and {Am^l ~ 2.45xlCT 3 cV 2 [1], 
where Am| 3 > or Am^ < refers to normal or in- 
verted mass hierarchy spectrum respectively. One of the 
most famous approaches to generate the active neutrino 
masses is the so-called "Type-I seesaw mechanism", in 
which one adds N right-handed neutrinos Nr ( (i = 1 — N) 
to the SM and the active neutrino masses can be obtained 
by block diagonalizing the mass matrix of left and right 
handed neutrinos, 

"^a = -L MRi ■ (1) 

Here a, = e, (jl,t represent the flavor indices of the SM 
fermions, Md is the Dirac mass matrix formed through 
the Yukawa interactions between left- and right-handed 
neutrinos, and Mr 4 are the Majorana masses of right- 
handed neutrinos. Since the right-handed neutrinos are 
completely neutral under the SM gauge symmetries, the 
Majorana mass Mr is a gauge invariant quantity and 
are often termed "sterile" neutrinos. At least two ster- 
ile neutrinos are needed to accommodate the two mass 
splits observed experimentally. The mass scales of Mo ai 
and are free parameters and cannot be fixed by os- 
cillation experiments alone. 

There are, however, results from the LSND [2] and 
the MiniBooNE [3] which cannot be accommodated in 
three active neutrinos description and may need to in- 
troduce one or more sterile neutrinos at the eV scale to 



* cmanshu® phy s . sinica . edu . tw 

* ryo. takahasi88@gmail . com 



fit the data (e.g. see [4, 5]). It should be mentioned that 
the light sterile neutrinos are also welcome in order to 
have successful Big Bang Nucleosynthesis (BBN) [6, 7], 
and are consistent with the preference for additional 
relativistic degrees of freedom, N e ff = 4.34^0.88' ob- 
served from the current cosmic microwave background 
(CMB) anisotropy probe and the large-scale structure 
(LSS) data [1]. 

Meanwhile, there is an increase in the amount and 
precision of cosmological data indicates that about 80% 
of the matter content in the universe is non-baryonic 
dark matter (DM) . The study of nature of DM is one of 
the main topics in cosmology, astrophysics, and particle 
physics. The Cold Dark Matter (CDM) is widely studied 
partly because the WIMP (weakly interacting massive 
particle) may reveal its signal at LHC and is predicted 
in many popular models (supersymmetric models, etc). 
However, it has been noticed that a sterile neutrino with 
mass at the keV scale and with small mixing to the active 
neutrinos can make up the DM in the form of Warm Dark 
Matter (WDM) [8, 9]. Additionally, it was pointed out 
that the keV sterile neutrino might play an important 
role in explaining the pulsar kicks [10]. 

Finally, the level of one out of ten billions excess in 
the amount of matter over antimatter is a long stand- 
ing puzzle for high energy physicists. The observation 
hints that baryon number (B) and/or lepton number (L) 
are violated in certain physical processes. Grand uni- 
fied theories (GUTs) naturally provide a framework for 
breaking B and L, in which the fundamental fermions 
- quarks and leptons - are arranged in the same multi- 
plets, and the out-of-equilibrium decays of heavy gauge 
bosons or colored Higgs bosons Hq will generate the suf- 
ficient baryon asymmetry around the scale of grand uni- 
fication [11-15]. It was then recognized that the Stan- 
dard Model (SM) violates B + L symmetry through the 
SU(2)l global anomaly [16]. The process is not sup- 
pressed during the period 100 GeV < T < 10 12 GeV, and 
the solution is called " sphalerons" [17]. The sphaleron ef- 
fect violates B + L but conserves B — L, and therefore, it 
would erase any primordial B + L asymmetry. We notice 
that any grand unification theories with higher symme- 
tries respects B — L symmetry. For example, B — L is a 
global symmetry for SU (5) GUTs and a local symmetry 
for 5O(10) GUTs respectively. The GUT-baryogenesis, 
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TABLE I: The content of three sterile neutrinos models 
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therefore, is not able to explain baryon asymmetry in our 
universe (BAU). To solve the problem one has to generate 
B — L asymmetry by violating pure baryon number [18] or 
by violating pure lepton number [19] (leptogenesis, e.g.), 
and the sphaleron process will convert partially B — L 
asymmetry into baryon asymmetry. We adopt the con- 
struction that one heavy sterile neutrino causes lepton 
asymmetry during the epoch the sphalerons are ineffec- 
tive, and the late decay of colored Higgs will generate the 
observed BAU. 

We consider a scenario of three sterile neutrinos 
Nn i(i=1 _ 3) with hierarchical mass spectrum (M Rl ~eV, 
Mr 2 -keV, M Rs ~ 10 13 ~ 15 GeV), in which the lightest 
one may help to explain the neutrino oscillation anoma- 
lies, keV-scale sterile neutrino is the candidate of dark 
matter, and the heaviest state Nr 3 would resurrect the 
GUT-baryogenesis. Three sterile neutrinos can be intro- 
duced to cancel the additional gauge anomaly for any 
theory beyond SM with extra gauge U(1)b-l symme- 
try. We show this hierarchical spectrum of sterile neuri- 
nos simultaneously satisfy the observations, and how our 
scenario fits in the framework of GUT theories. 

It has been proposed that models of SM with (three) 
additional sterile neutrinos are phenomenologically vi- 
able [20-22]. The so-called j/MSM (v Minimal Standard 
Model) [20], in which a mass of a keV sterile neutrino 
is responsible for DM, and two heavier states with de- 
generate masses lain in the range 1 GeV <~ 100 GeV are 
required to be in thermal equilibrium around electroweak 
scale in order to generate BAU through the resonant neu- 
trino oscillations. The split seesaw model with three ster- 
ile neutrinos living in the extra dimension (ED) is shown 
to be able to solve DM and BAU as well [21]. By uti- 
lizing an exponential factor in the size of ED one can 
split the Major ana masses of N Ri with relative mild pa- 
rameters associated to their locations in ED. Recently a 
flavor symmetry model [22] proposed by Barry, Rodejo- 
hann, and Zhang (BRZ), it consists of two N Rl 2 masses 
at eV-scale and one Nr 3 at keV-scale. The two eV-scale 
sterile neutrinos are used to explain LSND and Mini- 
BooNE anomalies while the keV sterile neutrino is the 
WDM particle. The scenarios are summarized in Ta- 
ble I. These setup can answer two of the three puz- 
zles we mentioned above while our hierarchically act- 
ing sterile neutrinos (HASN) scenario would explain the 
three puzzels simultaneously. The splittings of the sterile 
neutrino masses can be achieved by implementing split 
seesaw mechanism [21] or Froggatt-Nielsen (FN) mecha- 



nism [23] to the model. 

The Lagrangian which is relevant to neutrino masses 
has the form 



C = £ S m + iN Ri jBN Ri - y ai HH a N Ri 

Mr - 

7r-N c Ri N Ri +h.c. 



(2) 



Here £sm is the SM Lagrangian, l a are SU(2)l leptonic 
doublets with flavor index a, H is SM Higgs, y a i are the 
Yukawa couplings, and c is charged conjugation. The 
Majorana mass matrix of sterile neutrinos is chosen to be 
diagonal without loss of generality. The 6x6 neutrino 
mass matrix is given in the form 



M D 

Ml Mr 



(3) 



in the basis (v e ,v f _ L ,v T ,N Rl ,NR 2 ,NR 3 ), and Mr = 
diag(0(eV),0(keV),0(lO 13 - 15 )GeV). Here we give a 
brief comment on a realization of such a hierarchical mass 
spectrum of right-handed neutrinos. One of simple exam- 
ples to realize it is to utilize the split seesaw mechanism, 
in which the spinor fields are introduced in a flat five di- 
mensional (5D) spacetime whose compactification length 
of extra dimension is £ and all SM particle are assumed 
to live in a 4D-brane. After solving the 5D Dirac equa- 
tion and identifying the zero-modes of the 5D spinors 
with the right-handed neutrinos, the effective (4D) right- 
handed Majorana masses are described by exponential 
functions as M Ri = 2K, i miV B - L /(M(e 2mie ' — 1)) where 
Ki, nii, «b-l, and M are a coupling constant of order 
one, bulk masses for 5D spinors, U (1) b-l breaking scale, 
and 5D fundamental scale respectively. In this mecha- 
nism, one can easily obtain a hierarchical right-handed 
neutrino mass spectrum such as {Mr 1 , Mr 2 , Mr 3 ) 

I eV, 1 keV, 10 13 GeV) within a set of moderate param- 
eters when one takes k, = 1, «b-l = 10 15 GeV, and 
(Me,m 1 t,m 2 l,m 3 £) = (30,27.9,24.4,1.03) as reference 
values. The FN mechanism can also give a hierarchical 
mass spectrum with appropriate C/(1)fn charges. 

After electroweak symmetry breaking where Higgs de- 
velops its vacuum expectation value (VEV) v = 174 GeV, 
one gets Dirac neutrino mass terms. The left-handed 
neutrinos receive their Majorana masses through seesaw 
mechanism, we obtain 
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(4) 



m so i ^ 



tog 
M Rl 



Mr 2 



for both NH and IH, (5) 

(6) 



A/ 



for NH 
for IH 



at the leading order, where NH and IH mean the nor- 
mal hierarchy and inverted hierarchy respectively. The 
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indices a and ft in m„ 3 for NH should correspond to only 
fx and r in order to be consistent with the current data 
of neutrino oscillation experiments, that is, there are a 
maximal atmospheric, a large solar, and a small reactor 
mixing angles. By choosing a set of appropriate values of 
the Yukawa couplings, the experimentally observed mix- 
ing angles can be always fitted in our scenario. The de- 
generated mass spectrum of active neutrinos can be also 
realized. 

The dark matter candidate of the scenario is decay- 
ing DM. The keV sterile neutrino Nr 2 should live longer 
than the age of the universe and can be estimated as 
t Nr2 ~ 5 x 10 26 (MR 2 /kcV)- 5 (lCr 8 /e 2 )s, here 9 is the 
mixing between keV sterile neutrino and active neutri- 
nos. The generic way to produce DM is through the 
active-sterile neutrino oscillations [8], however, the abun- 
dance is constrained by the X-ray observations [24] (also 
see [25] and references therein), structure formation sim- 
ulations, and the Lyman-a bounds [26]. One way to re- 
lax the restrictions was proposed by Shi and Fuller [9] 
that an enhancement of the production of keV sterile 
neutrino can be realized via lepton-number-driven MSW 
(Mikhcycv-Smirnov-Wolfenstein) effect. The other pos- 
sibility is the Nr 2 pair production via C/(1)b-l gauge 
boson exchange [21]. It has been shown that as long as 
the reheating temperature is about 10 13 GeV one can 
account for the relic abundance of DM. The correspond- 
ing Yukawa couplings of sterile neutrino DM for the re- 
quired mass 0(1) keV< Mr 2 < O(10) keV are typically 
restricted to 0(1(T 15 ) < \y a2 \ < 0(1O~ 13 ) to satisfy 
astrophysical constraints (see e.g. [25] and references 
therein) . This means that terms from the sterile neutrino 
DM through the seesaw mechanism does not contribute 
to the atmospheric and solar neutrino mass scales shown 
in (4)- (6). The Yukawa couplings for the 1st and 3rd 
generations of right-handed neutrinos are approximated 
as_|2/* 3 2//33| 1/2 ~ 0(0.1) and \y* al y f n\ 1/2 ~ 0(10^- 
10 12 ) to satisfy the atmospheric and solar scales with 
(M Rl ,M R;s ) = (1 cV, 10 13 GeV) respectively It is seen 
that the construction of active neutrino mass spectrum 
in HASN scenario is consistent with the constraints on 
keV sterile neutrino DM. It can be also found that the 
atmospheric scale can be derived from the ratio of the 
right-handed neutrino mass, Mr 3 <~ 0(1O 13 ) GeV, and 
the corresponding Dirac masses, \ya3VP3\ 1 ^ 2y ~ O(10) 
GeV, when Nr 3 gets integrated out. While the solar scale 
comes from the seesaw relation between M Rl ~ 0(1) 
eV and \y* al ypi\ ll2 v ~ 0(0.1) eV. Finally, the ratio of 
mass scales between the 2nd generation of sterile neu- 
trino (DM), Mr 2 <~ 0(1) keV, and the corresponding 
Dirac masses, \y* a2 V r /32 1 ] L/ ' 2 « ~ 0(1O _3 -1O _1 ) eV, is too 
steep to contribute to the active neutrino mass (atmo- 
spheric and solar) scales. Therefore, the sterile neutrinos 
are hierarchically acting also for giving the active neu- 
trino mass scales. The Yukawa structure realizing the 
scenario can be obtained in both split seesaw and FN 
mechanisms with appropriate model parameters. 

Now we come to the phenomena of neutrino oscillation 



anomalies. The LSND — > D e transitions anomaly re- 
ported a 3.8cr excess of v e candidate events, in which the 
neutrino fluxes were produced by dumping 800 McV pro- 
tons into a "beam stop" which mostly generate 7r + , and 
neutrinos (anti-neutrinos) are the decay products of pi- 
ons. The probability that v a oscillates into v\, is given by 
P(ab) = sin 2 (20) sin 2 (1.27Am 2 -§), where 9 is the mixing 
angle, L is the neutrino travel distance in the unit of me- 
ter, and E is the neutrino energy in MeV. The typical 
anti-neutrinos energies are a few MeV for reactor exper- 
iments, the excess is interpreted as the hints for — > v e 
oscillation with Am 2 <~ 1 eV 2 . This indicates at least 
one sterile neutrino with mass at the eV-scale. Then the 
MiniBooNE experiment set out to check the excess events 
in the — > v e transitions and found the parameters were 
not compatible with LSND [27]. However, more recently 
the MiniBooNE accumulated more anti-neutrino oscilla- 
tion data and reported the excess electron anti-neutrino 
appearance is reconciled with LSND results [3]. To ac- 
commodate the neutrino and anti-neutrino data the ad- 
ditional CP violation has to be invoked. One simple way 
is to add two sterile neutrinos at eV scale (the so-called 
(3+2) scheme) to neutrino sector, the CP violation at 
short-baselines would let to reconcile both LSND and 
MiniBooNE results [4, 5, 28-30]. In our scenario a (3+1) 
scheme together with nonstandard interactions (NSI) of 
neutrinos will allow to fit the data [4] . The new interac- 
tions may modify the charged and neutral currents, and 
provide the new sources of CP violation. They may affect 
the neutrino oscillations via the production, propagation, 
and the detection processes. The four-fermion operators 
can be expressed at low energies as 

£nsi - ZV2G F J2 ^y ^Lal^LfsKh^fL^) + h.c.,(7) 

/ 

where Gf is Fermi constant, / represents fermions 
(charged leptons and quarks), and a, ft are flavor indices, 
and L, R are chiralities. The new interactions can be in- 
duced from several possibilities of physics beyond SM. 
For example, in GUT theories the ip(10) © ip(5) fermion 
representations of SU(5) are coupled to a 5(Hc) and a 
5(Hc) representation of Higgs, and the Yukawa interac- 
tions rcad£ y = ip(10) T \ u ip(10)H c +ip(10) T \ d ip(h)H c . 
It has been shown that one can fit to global short-baseline 
data for e afs ~ 0(1O~ 2 ) [4]. 

In the context of our consideration, we adopt that 
a sterile neutrino Nr 3 is heavier than the colored 
Higgs bosons, and a lcpton number violating interac- 
tion lcj>lcj>/(2Mii 3 ), here <f> is the Higgs doublet, keeps 
in thermal equilibrium when T > 10 12 GeV. There- 
fore, all lepton asymmetry generated by colored Higgs 
decay is erased by the process I + <f> — > I + <ft while 
the generated baryon asymmetry remains intact, and so 
B — L =/= can be satisfied. When temperature drops be- 
low 10 12 GeV the sphaleron transitions become effective, 
as the results, the produced baryon asymmetry is par- 
tially converted into the lepton asymmetry but a resid- 
ual baryon asymmetry remains, and thus the observed 
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BAU can be generated. In the case of SU(5) GUT wc 
have the Yukawa couplings ip(10) T \ u( - k 1 ip(W)H^ ] and 
ip{lG) T X d ^^{b)H^ ] (k = 1,2). The size of BAU is 
calculated as [15, 19, 31, 32] 

Ja , s M_a.„, s ,,, r! ._i_, (!| 

where K s I£l T _ m „ ~ GcV f -Y' 2 is the 

washout factor and e B ~ ^ • 10" 2 [F(a;) - F(l/a;) + 
G(a:) — G(l/a;)] is the CP-asymmetry with T, iJ, mjj c , 
g*, ir are the decay rate, expansion rate, colored Higgs 
mass, degrees of freedom g*\T~m Hc — 53, mass ratio 
m 2 (2) /m 2 (1) , respectively. The functions F and G are 

defined as F(x) ~ 1 - a; In (±±*) and G{x) = -^j. The 
factor 0.35 comes from the sphaleron process, and we take 
771 = sin(arg[tr(A d ( 1 )tA d ( 2 )A"( 1 )tA u ( 2 ))]) and m ~ r/ 2 . It 
is seen that we can realize Yab — 8.75 x 10~ n when 
we set (m H (i),m H (2)) = (9 x 10 12 ,8 x 10 12 ) GeV and 



rji ~ —0.444. These values are consistent with this baryo- 
genesis scenario and above discussion of the right-handed 

neutrinos mass spectrum, that is, 10 12 GeV < m <;> < 

n c 

Mr 3 < 10 15 GeV. A heavier mass of the corresponding 
sterile neutrino as 10 14-15 GeV is also possible for this 
baryogenesis. 

The hierarchical spectrum of sterile neutrinos (eV, 
keV, 10 13 ~ 15 GeV) is a simple and economical scenario, 
especially it can be embedded in many frameworks be- 
yond SM. In light of the puzzles from neutrino oscilla- 
tion anomalies, dark matter, and baryon asymmetry of 
the universe (BAU) we have shown that this scenario is 
phenomcnologically viable. Searching for deviations from 
standard three active neutrino oscillations and the X-ray 
astronomy will offer opportunity to test the scenario. 
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